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Abstract.  In  this  paper  we  discuss  a  problem  of  generalization  of  binomial  distributed  triangle, 
that  is  sequence  A287326  in  OEIS.  The  main  property  of  A287326  that  it  returns  a  perfect  cube 
n  as  sum  of  n-tli  row  terms  over  k,  such  that  0  <  k  <  n  —  1  or  1  <  k  <  n,  by  means  of  its 
symmetry.  In  this  paper  we  have  derived  a  similar  triangles  in  order  to  receive  powers  m  =  5,  7  as 
row  items  sum  and  generalized  obtained  results  in  order  to  receive  every  odd-powered  monomial 
n2m+i,  m  >  0  as  sum  of  row  terms  of  corresponding  triangle.  In  other  words,  in  this  manuscript 
are  found  and  discussed  the  polynomials  Dm(n,k)  and  Um{n,k),  such  that,  when  being  summed 
up  over  k  in  some  range  with  respect  to  m  and  n  returns  the  monomial  n2m+1. 
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1.  Structure  of  the  manuscript 


The  problem  of  finding  expansions  of  monomials,  binomials,  trinomials,  etc.  is  classical  and  a  lot 
of  theorems  have  been  found,  the  most  prominent  examples  are  Binomial  Theorem  [2],  Multinomial 
theorem,  Worpitzky  Identity  [30],  [31],  Stirling  numbers  of  second  kind  and  falling  factorial  identity, 
[36 j .  etc.  In  this  paper  we  try  to  solve  the  classical  problem  of  finding  expansions  of  monomials. 
We  start  from  binomial  distributed  triangle  A287326,  HU  in  OEIS.  The  main  property  of  A287326 
that  it  returns  a  perfect  cube  n  as  re-th  row  sum,  starting  from  0, ...,  n  —  1  or  from  1, ...,  n  by  means 
of  its  symmetry.  Therefore,  the  following  question  stated: 


•  Is  there  a  generalization  of  A287326  in  order  to  receive  monomial  nt,  t  >  3  as  sum  of  n-th 
row  terms  of  corresponding  triangle,  where  t  is  natural? 

Finding  an  analogs  for  t  =  5,7  in  section  [3]  we  answer  to  above  questions  positively.  Could  this 
process  be  continued  for  each  odd  t  =  1, 3,  5,  7...  similarly?  Positive  answer  to  this  question  is  given 


by  theorem  (3.30) 
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2.  Introduction 


Let  describe  the  derivation  of  the  sequence  A287326  in  OEIS.  Sequence  A287326  returns  the 
perfect  cube  n  as  row  sum  over  k,  0  <  k  <  n  —  1,  as  well  as  sum  over  1  <  k  <  n,  by  means  of  its 
symmetry.  First,  consider  a  difference  table  of  perfect  cubes  (0,  eq.  7) 


(2.1) 


n 

A u(n3) 

Ai(n3) 

A  z(nA) 

Aa(n3) 

0 

0 

1 

6 

6 

1 

1 

7 

12 

6 

2 

8 

19 

18 

6 

3 

27 

37 

24 

6 

4 

64 

61 

30 

6 

5 

125 

91 

36 

6 

6 

216 

127 

42 

6 

7 

343 

169 

48 

6 

8 

512 

217 

54 

9 

729 

271 

10 

1000 

Table 


1:  Difference  table  of  perfect  cubes  n,  0  <  n  <  10  up  to  3rd  order. 


Reviewing  above  table,  we  have  noticed  that 


(2.2) 


A(03) 

=  1  +  6  ■ 

'  »  =  6(f)  +  ffl 

A(l3) 

=  1  +  6' 

0  +  6- 1  =  6©  +  © 

A(23) 

=  1  +  6  ' 

0  +  6-1  +  6-  2  —  6(2)  +  (q) 

A(33) 

=  1  +  6' 

■0  +  6-1  +  6-  2  +  6-  3  —  6(2)  +  (!j 

A(n3)  =  l  +  6-0  +  6-l  +  6-2  +  ---  +  6-n  =  6(n+1)  +  (n+1) 

Above  difference  identity  is  closely  related  to  Faulhaber’s  sum  of  cubes,  where  n3  =  6("i^1)  + 
see  ([2U,  P-  9)-  Note  that  A 2(ra3)  could  be  found  similarly  using  above  identity,  i.e  A2(n3)  = 

6(S)  +  M- 

Property  2.3.  (Generalized  finite  difference  of  power  using  Faulhaber’s  formula).  Consider  the 
identities,  m,  p.  9).  For  every  odd  power 


f"1  =  (?) 

n3  =  6(^)  +  (?) 

n5  =  120  (n+2)  +  30(n+1)  +  (?) 


n2m-i  =  ^  (2fc-l)!T(2m,2fc)(n+fcT11) 

<  1  <k<m 

The  coefficients  in  these  formulas  are  related  to  what  Riordan  [22]  has  called  central  factorial  num¬ 
bers  of  the  second  kind.  In  his  notation, 

xm=  T(m,k)xW,  zM  =x(x  +  |-l)(®  +  |-2)--.(x+|  +  l) 

1  <k<m 

The  coefficients  T(2m,2k)  are  always  integers,  because  the  x^k+2^  =  x^(x2/k4)  implies  the  recur¬ 
rence 

T(2m  +  2,  2k)  =  k2T(2m,  2k)  +  T(2m,  2k -2). 


SERIES  REPRESENTATION  OF  POWER  FUNCTION 


3 


We  can  find  the  first  order  finite  difference  of  odd  power  as  decreasing  the  variable  of  corresponding 
binomial  coefficients  by  1,  for  example 

An1  =  ffl 

A"3  =  6(”+‘)  +  ffi 

An5  =  120("J'2)  +  30("J‘)  +  (;) 


A?r 


2m—  1 


=  E  (2k-l)\T(2m,2k)(n+k_-21) 

l<k<m 


Continue  similarly,  we  can  express  each  difference  of  order  t  >  1.  The  central  factorial  numbers 
of  the  second  kind  (2 k  —  l)!T(2m,  2k)  in  above  identities  are  terms  of  OEIS  sequence  A303675  and 
generated  by 


(2.4) 


(2k  -  l)\T(2m,  2k)  =  -  £(-1^  (r  -  j?n  =:  Vm,k, 
r  j= o  3  J  ' 


where  r  =  n  —  k  +  l,  this  formula  was  provided  by  Peter  Luschny  in  m-  Repeated  sums  are  equally 
easy,  in  Knuth’s  notation  (see  [21],  p.  10) 

n  +  k  +  r 


r_2m+l 


Tin 


^  )  Vm, k 


2k  —  1  +  r 


1  <k<m 

Therefore,  reviewing  the  difference  as  inverse  operator  to  summation  for  every  odd  t  >  0  and  m  >  0, 
we  have  identity 

2m+l  _  V-  „  fn  +  k-r' 


Arn 


T  * 


m.k 


1  <k<m 


2k  —  1  —  r 


□ 


By  property  2.3  we  rewrite  the  cubes  as 
(2.5) 


n3  = 


£  6lkt1)  +  rk 

l<k<n 


Rewrite  above  expression  with  set  every  binomial  coefficient  to  be  (r£)  =  1  +  2  +  •  •  •  +  (n  +  1), 
then 

n3  =  (1  +  6  •  0)  +  (1  +  6  •  0  +  6  •  1)  +  •  •  •  +  (1  +  6  •  0  +  •  •  •  +  6  •  (n  -  1)) 
Particularizing  above  expression,  we  get 

(2.6)  n3  =  n  +  (n  —  0)  •  6  •  0  +  (n  —  1)  •  6  •  1  +  •  •  •  +  (n  —  (n  —  1))  •  6  •  (n  —  1) 


Provided  that  n  is  natural.  Now,  let  apply  a  compact  sigma  notation  on  (2.6),  thus 

+  £  6 k(n  —  k ) 


(2.7) 


■? 

n  =  n 


1  <k<n 


As  sum  Ei <k<n^k(n  —  k )  consists  of  n  terms,  we  have  right  to  move  n  in  (2.7)  under  sigma 
notation,  we  get 

(2.8) 


n3  = 


£  6 k(n  —  k)  +  1 


l<fc<n 
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Property  2.9.  (Proof  of  symmetry).  Let  be  a  sets  A(ri)  :=  {1,  2,  n},  £?(n)  :=  {0,  1,  n}, 

C(n )  :=  {0,  1,  . . . ,  n  —  1},  Zet  6e  expression  $2. tty  denoted  as 

M(n,  C{n))  =  ^  6k(n  —  k)  +  1 

fceC(n) 

where  n  is  natural-valued  variable  and  C(n)  is  iteration  set  of  then  we  have  equality 

(2.10)  M(n,  A(n ))  =  M(n,  C(n )) 

Let  review  and  denote  expression  \2.6 ])  as 

U(n,  C(n ))  =  n  +  6-  ^  k(n  —  k) 

kdC{n) 


then 


(2.11)  t/(n,  A(n))  =  U(n,  B(n))  =  U(n,  C(n )) 

Other  words,  changing  of  iteration  sets  of '{2. 6])  and  (f£.<S|)  by  A(n),  B(n),  C(n )  and  A(n ),  C{n), 
respectively,  doesn't  change  resulting  value  for  each  natural  x. 

Proof.  Let  be  a  plot  y{n ,  fc)  =  6£;(n  —  fe)  +  1,  k  6  M,  0  <  k  <  10,  given  n  =  10 


Figure  2.  Plot  of  6fc(n  —  fc)  +  1,  k  €  M,  0  <  k  <  n,  where  n  =  10. 


Obviously,  being  a  parabolic  function,  it’s  symmetrical  over  ” ,  hence  equivalent  M(n,  A(n))  = 
M(n,  C(n))  follows.  Reviewing  (2.6)  and  denote  u(n,k )  =  kn  —  k 2,  we  can  conclude,  that  u(n, 0)  = 
u(n,n )  =  0,  then  equality  of  t/(n,  .A(n))  =  U(n,  B{n ))  =  17(n,  C{n))  immediately  follows.  This 
completes  the  proof.  □ 


Review  above  property  (2.9).  Let  be  an  example  of  triangle  built  using 


Definition  2.12.  For  every  n  >  0 

(2.13)  Di(n,  k)  =f  6 k(n  —  k)  +  1,  0  <  k  <  n 


over  n  from  0  to  n  =  4,  where  n  denotes  corresponding  row  and  k  shows  the  item  of  row  n. 
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(2.14) 


Row  0 
Row  1 
Row  2 
Row  3 
Row  4 


1 

1  1 
1  7  1 

1  13  13  1 

1  19  25  19  1 


Figure  3.  Triangle  generated  by  Di(n,k ),  sequence  A287326  in  OEIS,  [TTj. 


Note  that  n-th  row  sum  of  Triangle  (2.14)  over  0  <  k  <  n  —  1  returns  perfect  cube  n.  We  can  see 


that  each  row  with  respect  to  variable  n  =  0,  1,  2,  3,  4, ...,  has  Binomial  distribution  of  row  terms. 


One  could  compare  Triangle  (2.14)  with  Pascal’s  triangle  in,  m. 


Row  0 
Row  1 
Row  2 
Row  3 
Row  4 


1 

1  1 

1  2  1 
13  3  1 

1  4  6  4  1 


Figure  4.  Pascal’s  triangle,  sequence  A007318  in  OEIS,  [lj. 


Let  us  approach  to  show  a  few  properties  of  triangle  (2.14)  and  L\(n,k). 


Properties  2.15.  Properties  of  triangle  (2.14) 


(1)  Summation  of  n-th  row  of  triangle  (2.14)  over  k  from  0  to  n—  1  returns  perfect  cube  n  >  0 
as  follows 


(2.16) 


to3  = 


Y  Di(n,k) 


1  <k<n 


(2)  First  item  of  each  row’s  number  corresponding  to  central  polygonal  numbers  sequence  a(n)  = 
n  +2+2  (sequence  A000124  in  OEIS,  [TSjj  returns  finite  difference  of  consequent  perfect 


cubes.  For  example,  let  be  a  k-th  row  of  triangle  (2.14),  such  that  k  =  n  +2n+  ,  n  =  0,1,  2, ..., 


then  item 


(2.17) 


A(n3 4 5)  =  Di 


rr  +  n  +  2 


(2.18) 


(2.19) 


(3)  Items  of  (2.14)  have  Binomial  distribution  of  rows. 

(4)  Linear  recurrence,  for  every  k  and  n  >  0 

2D\  (to,  k)  =  D\ (to  +  1,  k)  +  D\  (to  —  1,  k) 

This  linear  recurrence  is  direct  result  of  second  order  binomial  transform  of  D\(n,  k)  by  n. 

(5)  Linear  recurrence,  for  each  n  >  k 

2D\  (to,  k)  =  D\  (2 n  —  k,  k)  +  D\  (2 to  —  k,  0) 


(6)  From  ( 1.24)  for  every  n  >  0  follows 

(2.20) 


TO3  = 


£  Di {n,k)=  Y.  Ap2  +  t  +  2.l 

Kk<n  l<k<n 


(7)  Triangle  (2.14)  is  symmetric,  i.e 

(2.21) 


D\  (to,  k)  =  D\  (to,  to  —  k) 
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Corollary  2.22.  Review  identity  (2.16)  in  sense  of  summation  of  D\(n,k)  over  some  range  W 
with  max(W)  =  T,  then  (2.16 ])  returns 


Y  D1(n,k)  =  U1(T,0)n°-U1(T,l)n1, 

1  <k<T 


where  T  =  1,2,3,...,N.  By  property  we  rewrite  above  expression  as 


Y  Di(n>  k )  =  Ui(T-  l,0)n°  -  Ui  (T  -  1,  l)n1 

0<k<T—l 


Below  we  show  a  few  initial  terms  of  U\(T  —  1,  j),  Ui(T,j),  j  =  0, 1  coefficients 


T 

U\(T  —  1, 0) 

Ui(T-  1,1) 

Ui(T,  0) 

Ui(T,  1) 

1 

1 

0 

-5 

6 

2 

-4 

6 

-28 

18 

3 

-27 

18 

-81 

36 

4 

-80 

36 

-176 

60 

5 

-175 

60 

-325 

90 

6 

-324 

90 

-540 

126 

7 

-539 

126 

-833 

168 

8 

-832 

168 

-1216 

216 

9 

-1215 

216 

-1701 

270 

10 

-1700 

270 

-2300 

330 

Table  5.  Table  of  coefficients  U\{T  —  1,  0),  U\(T  —  1, 1),  U\{T ,  0),  Ui(T,  1)  given  T  =  1, 10. 
Therefore,  for  every  integer  n  >  1,  T  =  n, 

' U\{T,  0)n°  —  U\{T,  l)nl , 

Ui(T-  l,0)n°  -  U\(T  —  1,  l)n1. 

Coefficients  |Hi(T,  0)|  are 


VT  =  n  :  rr  = 


Coefficients  \U\{T  —  1, 1)|  are  terms  of  sequence  A028896  in  OEIS, 
terms  of  sequence  A275709  in  OEIS,  [20] . 


□ 


In  this  section  we  have  derived  a  binomial  distributed  triangle  (2.14),  such  that  perfect  cube 
n  could  be  found  as  sum  of  n-th  row  terms  of  triangle  (2.14)  over  k,  in  ranges  1  <  k  <  n  or 
0  <  k  <  n  —  1,  where  n  is  natural.  Therefore,  the  follow  question  is  stated: 

Question  2.23.  Is  there  a  generalization  of  A287326  in  order  to  receive  monomial  nt,  t  >  3  as 
sum  of  n-th  row  terms  of  corresponding  triangle,  where  t  is  natural? 

3.  Generalization  of  sequence  A287326 


In  order  to  get  analogs  of  Triangle  (2.14)  one  should  solve  a  system  of  equations,  where  unknowns 
are  coefficients  of  polynomial  and  variable  of  polynomial  is  k(n  —  k).  Triangle  (2.14)  is  generated 
by  polynomial  D±(n,  k ),  n  >  0,  0  <  k  <  n,  defined  by  (2.12).  Here,  let  derive  a  triangle  generated 
by  polynomial  D-2{n,k),  n  >  0,  0  <  k  <  n,  such  that  sum  of  n-th  row  terms  over  k,  in  ranges 
I<fc<nor0<fc<n-1  returns  n5,  where  n  is  natural 

Example  3.1.  We  suspect  that  n-th  row  of  triangle  that  returns  n5  as  sum  of  n-th  row  terms  over 
A;,  in  1  <  k  <  n  or  0  <  k  <  n  —  1  is  generated  by 

(3.2)  T>2(n,  k )  =  A2t2{n  ~  k)2k2  +  H.2;i(n  —  A;)1^1  +  A2${n  —  k)°k°,  n  >  0,  0  <  k  <  n, 
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where  A2t2,  A-2yi,  A2fi  are  unknown  coefficients.  Assume  that  for  every  integer  n  >  0  holds 
(3.3)  n5  =  22  D2(n,k). 


1  <k<n 


To  determine  the  coefficients  ^4.2,2,  ^4.2,1,  ^4.2,0?  in  (|3.2|)  let  rewrite  (|3.3[)  in  extended  view 
(3.4)  A2,2  22  k2(n  -  k)2  +  A2,i  k(n  -  k)  +  22  A2)o 

l<fe<n  1  <k<n  1  <k<n 

=  A2  2  22  k2(n2  -  2nk  +  k2)  +  A2ii  22  kn  -  k2  +  22  A2,o 
l<fc<n  1  <k<n  1  <k<n 

=  A2j2  22  k2n2  —  2nk3  +  kA  +  A2.i  22  k n  ~  k2  +  22  ^2,0 

l<fc<n  l<k<n  1  <k<n 

=  A2^2n2  22  k2  -  2A2j2n  22  k3  +  A2j2  22  k4  +  A2,in  22  k 

1 <k<n  1 <k<n  1 <k<n  1 <k<n 

—  a2)i  '*22  k2  +  22  ^2>° =  n5- 

1 <k<n  1 <k<n 

Thus,  we  have  received  expression  containing  sums  of  powers  of  successive  natural  numbers,  where 
powers  are  {1,  2,3,4}.  These  formulas  contain  so-called  Bernoulli  numbers,  |14i.  For  mentioned 
powers  of  successive  natural  numbers  formulas  are  following 


(3.5) 

(3.6) 

(3.7) 

(3.8) 


E  * 

l<fc<n 

E  * 

l<fc<n 

E  * 

l<fc<n 

E  * 

l<fc<n 


n2  +  n 


2  2n3  +  3n2  +  n 
~  6  ’ 

3  n4  +  2  n3  +  n2 


4  6n5  +  15n4  +  10n3  —  n 

30  ' 


Now  we  substitute  above  identities  to  (|3.4|),  respectively, 
A2y2n 


2  2n3  +  3  n2  +  n 


6 


,  n4  +  2  n3  +  n2  .  6n5  +  15n4  +  10n3  —  n 

-  2 A2y2n - - - T  A2  2j - 


30 


+  A2\n 


n2  +  n 


.  2  n3  +  3n2  +  n 
—  A2j i - - - h  A2}0n 


2  ’  6 

Particularizing  the  elements  of  above  expression  and  moving  them  under  the  common  divisor,  we 
get 


(3.9) 


A2j2n5  —  A2y2n  +  30A2iq  .  n”  —  n 


30 


+  A: 


2,1  - 


6 


We  have  to  remember  that  expression  (3.9)  is  the  left  side  of  the  input  equation  (3.3).  Therefore, 
(3.10) 


5  A2  2n5  -  A2  2n  +  30A2,0  .  n3  -  n 

n  =  - - — b  A2  i - 

on 


30 
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In  order  to  satisfy  (3.10)  for  each  natural  n,  coefficients  ^2,0: ^2,1 
following  system  of  equations 


A2t2  should  be  a  solutions  of 


34^2,2 

=  1 

^4-2,1 

=  1 

30^2,0  —  A2l2 

=  0 

The  only  solution  of  above  system  is  ^2,2  =  30,  ^2,1  =  0,  =  1.  Hereby,  polynomial  D2(n,k ) 

takes  the  form 


(3.11)  T>2(n,  k)  =  A2y2(n  —  k)2k2  +  A2}i(n  —  k)3k3  +  A2fi(n  —  k)°k°  =  30 k2(n  —  k)2  +  1 
And  for  every  natural  n  >  1  holds 

(3.12)  n5  =  ^  D2(n,k)  =  ^  30k2(n  —  k)2  + 1 

1 <k<n  1 <k<n 


E 

0<fc<n— 1 


D2(n,  k)  = 


E 

0<A:<n— 1 


30k2  (n  —  k)2  +  1 


Let  show  few  initial  rows  of  triangle  built  by  D2(n,k ),  that  is  analog  of  triangle  (2.14),  which 
returns  monomial  n5  as  sum  of  n-th  row  terms  over  k,  asl<k<noi0<k<n  —  1,  by  means  of 
its  symmetry 


(3.13) 


1 

1  1 


1 

31 

1 

1 

121 

121 

1 

1 

271 

481 

271 

1 

1  480 

1081 

1081 

481 

1 

Figure  6.  Triangle  generated  by  polynomial  D2(n,  k),  n  >  0,  0  <  k  <  n,  sequence  A300656  in 

OEK,  H5]. 

Similarly,  finding  the  coefficients  A3io,  A3ji,  A3j2,  A3j3  in 

(3.14)  Z)3(n,  k )  =  A3)3k3(n  -  k )3  +  A3)2k2{n  -  k)2  +  A^i/c^n  -  A:)1  +  A3fik°(n  -  k)° , 
we  get  A3]3  =  140,  A3;2  =  —14,  A3j  =  0,  A3io  =  1,  therefore,  for  each  integer  n  >  0  holds 

(3.15)  n7  =  ^  D3(n,k)  =  ^  140k3 (n  —  k)3  —  14k2 (n  —  k)2  +  1 

l<k<n  1  <k<n 

=  ^2  D3(n,k)  =  22  140k3  (n  —  k)3  —  14k2  (n  —  k)2  +  1 

0<k<n-l  0<k<n—l 


Below  we  show  a  few  initial  rows  of  triangle  generated  by  polynomial  D3(n,  k ),  n  >  0,  0  <  k  <  n, 
the  analog  of  triangle  (2.14),  such  that  monomial  n7  could  be  found  as  sum  of  n-th  row  terms  over 
k,  as  1  <  k  <  n  or  0  <  k  <  n  —  1,  by  means  of  its  symmetry 
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(3.16) 


1 

1  1 


1 

127 

1 

1 

1093 

1093 

1 

1 

3793 

8905 

3793 

1 

1  8905 

30157 

30157 

8905  1 

Figure  7.  Triangle  generated  by  polynomial  1)3(71,  k),  n  >  0,  0  <  k  <  n,  sequence  A300785  in 

OEIS,  [T6]. 

We  assume  now  that  generalization  of  A287326  holds  for  odd  powers  of  the  form  2m  +  1,  m  = 
0, 1,  2, where  A287326  is  partial  case  for  m  =  1.  To  generalize  our  sequences  A287326,  A300656, 
A300785  for  every  odd  power  2m+l,  m  >  0  one  has  to  review  the  generating  functions  of  sequences 
A287326,  A300656,  A300785  as  follows.  Let  be  definition 


Definition  3.17. 

Dm(n,  k)  :=  Am,mkm(n  -  k)m  +  Am^m_1km~1  {n  -  -+ - b  Amfik° (n  -  k)° 

=  ^2  Am}jkj  (n  -  k)j , 

0<j<m 

where  Amj ,  0  <  j  <  m  are  unknown  coefficients. 


And,  we  assume  that 


(3.18) 


n2m+1  _  ^2  Dm(n,k). 

1  <k<n 


We  want  to  notice  that  as  we  used  a  compact  sigma  notation  on  definition  (3.17),  i.e  we  rewrite 
(3T7)  as  X)o <j<m.Am.,j  kJ(n  —  k)] ,  thus  the  sum  (13. 18h  returns  indeterminate  form  of  Dm(n,k )  = 
So <j<mJ^m,jk^{n  —  ky  on  step  k  =  n  as  Amfl(n  —  n)°k°  contains  the  term  (n  —  n)°  =  0°.  Some 
textbooks  leave  the  quantity  0°  undefined,  because  the  functions  x°  and  (P  have  different  limiting 
values  when  x  decreases  to  0.  But  this  is  a  mistake.  We  must  define 

Vx  :  x°  =  1, 

if  the  binomial  theorem  is  to  be  valid  when  x  =  0,  y  =  0,  and/or  x  =  y.  The  binomial  theorem 
is  too  important  to  be  arbitrarily  restricted!  By  contrast,  the  function  CP  is  quite  unimportant, 

For  example,  generating  functions 


for  A287326 
for  A300656 
for  A300785 


Note  that  Dm(n,k )  is  generalization  of  (2.12)  and  (3.11). 
of  sequences  A287326,  A300656,  A300785  are,  respectively 

D\(n ,  k)  =  1  +  6 k(n  —  k ), 

£>2(71,  k)  =  1  —  0 k(n  —  k)  +  30 k2{n  —  k)2, 

Ds(n,  k)  =  1  —  14 k(n  —  k )  +  0 k2{n  —  k)2  +  140fc3(n  —  k)3, 


Where  coefficients  Amj,  for  m  =  1,  2,  3  are  {Aij}j=0  =  {1,  6},  {A2j}2=0  =  {1,  0,  30},  {A3,j}3=0  = 
{1,  —14,0, 140}  in  definitions  of  generating  functions  of  A287326,  A300656,  A300785,  To  generalize 
above  result  in  order  to  receive  monomial  n2m+1  as  J21<k<n  Dm(n ,  k)  =  n2m+1,  m  =  0, 1,  2, ...  one 
has  to  solve  the  system  of  equations.  Complete  set  of  coefficients  {Am> 0,  . . . ,  Am;m}  such  that 
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Yli<k<n  Dm(n,  k)  =  n2m+1,  m  >  0  holds  can  be  found  by  solving  the  following  system  of  equations 

fDm{  1,0)  =  12™+1 
Dm(2,  0)  +  Dm(2, 1)  =  22m+1 

(3.19)  Dm(3,  0)  +  Dm{ 3, 1)  +  Dm{ 3,  2)  =  32m+1 

„Dm{r,  0)  +  Dm(r,  1)  H - b  Dm(r,r  -  1)  =  r2m+1,  r  >  m 

List  of  solution^]  of  system  (3.19)  is  split  and  assigned  to  OEIS  under  the  numbers  A302971 


(numerators  of  Amj)  and  A304042  (denominators  of  Amj).  To  reach  recurrent  formula  of  Amj, 
first  let  fix  the  unused  values  Amj  =  0,  for  j  <  0  or  j  >  m,  so  we  don’t  need  to  care  about  the 
summation  range  for  j,  then  by  expanding  (n  —  k and  using  Faulhaber’s  formula  [7],  we  get 

n— 1  n— 1  oo  /  .\ 

f3 


(3.20) 


—  kyw  =  EE 


k= o 

OO 


k= 0  i 


v  r  (~1}* 

\  %  J  i  +  j  +  \ 


l 

OO 


f ;(<+(+1W«+i-‘-b, 


i+j+l 


E(j)jyjyi(’  +  t  +  1)  ~  E  (() 


i,t 


t  ) 


o-i 


lj  1+  J  + 


(*) 


where  Bt  are  Bernoulli  numbers  p3j.  Now,  we  notice  that 

J  (2j+i; 


(3.21) 


W 


if  t  =  0; 


V  (  A  A  +  •?  +  A  =  J  (2j+i)(2/)  ’ 

^Lh+i  +  A  t  )  lfc^(  ift>„ 


In  particular,  the  last  sum  is  zero  for  0  <  t  <  j.  Now  we  revise  the  (*)  part  of  (3.20)  according  to 
results  of  (|3.21[),  thus 

1 


e(oa^iCtv+,-‘ 


(2j  +  1)(?) 


+  eA^».  i .  ,  wj+i-‘ 

t>  o 


2j  —  t  +  1 


Therefore,  (3.20[)  takes  the  form 


n—1 


(3.22) 


y~](n  —  kyW 


k= o 


(2j  +  1)(2/) 


+  E 


(-i)j 


t>0 


J  )  5#n2j+1_t 
2j-t  +  i; 


(*) 


AA  A1)*  p  j-i 
El*)*..  -  .  T^+i+i71 


i/i  +  i  +  l 


(o) 


Now,  we  keep  our  attention  to  (3.22)  and  we  have  to  remember  that  if  the  sum  over  some  variable  i 


contains  (•?) ,  then  instead  of  limiting  its  summation  range  to  i  =  0, ...,  j,  we  can  let  i  =  — oo, ...,  +oo 


Nolutions_system_2_4.txt  -  Mathematica  code,  provides  the  list  of  solutions  of  system  (3.191  up  to  r  =  11,  l24l. 
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since  (/)  =  0  for  i  outside  the  range  i  =  0,  ...,j  (i.e.,  when  i  <  0  or  i  >  j).  It’s  much  easier  to  review 
such  sum  as  summing  from  — oo  to  +oo  (unless  specified  otherwise),  where  only  a  finite  number 
of  terms  are  nonzero,  this  fact  is  discussed  in  [28]  as  well.  To  combine  or  cancel  identical  terms 


across  the  two  sums  in  (3.22)  more  easily,  we  introduce  t  =  2j  +  1  —  t  to  (*)  and  £  =  j  —  i  to  (o), 
respectively,  we  get 

77,-1 


(3.23) 


—  kyk> 

k= o 

OO 

£ 


1 


i=—oo 


(2j  +  l)(2/) 

j\  (-ir£ 


n 


2.7+1 


2j  +  1  -  l 


2j+l 


B2j+i-in 


pi+pff)"  +2„£ 


(-it 

2j  +  1  -  t  \l 


B2j+i-en 


Now,  using  the  definition  of  Am  j,  we  obtain  the  following  identity  for  polynomials  in  n 


(3.24) 


1 


(2j  +  1)  (2/) 


n 


27+1 


+  2  ^  Amj 


j,  odd  £ 


Rif  D  „( 


2  j  +  1  -  l 


=  n2m+1. 


Taking  the  coefficient  of  n2m+1  in  above  expression,  we  get  Amtin  =  (2m  +  1)(2™),  and  taking  the 
coefficient  of  x2d+1  for  an  integer  d  in  the  range  m/2  <  d  <  m  we  get  Am jd  =  0.  Taking  the 
coefficient  of  n2d+1  in  (3.24)  for  m/4  <  d  <  m/2  ,  we  get 

1  ^(2 m\  /  m  \  (—1)’ 

m 


(3.25) 


A 


m,d 


+  2(2  m  +  1) 


i.e 

(3.26) 


(2  d+l)ffl 

Am,d  =  (-  l)m_1 


(2m  +  1) 


2d  +  1/  2m  —  2d 

1 


B2m—2d  —  0, 


,  Bo.m  — 


2m— 2d- 


d\d\m\(m  —  2d  —  1  )\m  —  d 

Continue  similarly,  we  can  express  Amj  for  each  integer  j  in  range  m/2s+1  <  j  <  m/2s  (iterating 
consecutively  s  =  1,2,...)  via  previously  determined  values  of  Tm,d,  d  <  j  as  follows 


(-1) 


d- 1 


—  B2d-2j- 


(3.27)  +>3-<+  +  1)(/)fc£i'4^(2j/1) 

The  same  formula  holds  also  for  m  =  0.  Note  that  in  above  sum  m  have  to  be  m  >  2j  + 1  to  return 
nonzero  term  Amj. 

Definition  3.28.  We  define  here  a  generating  function  of  sequence  of  coefficients  Am  j,  such  that 


Ylk= o  Dm(n,  k)  — 


_  „2m+l 


n 


n  >  0,  m  >  0,  where  Dm(n ,  &:)  is  defined  by  (3.17) 


■ —  * 


0, 


if  j  <  0  or  j  >  m 


(2j  +  1)(  /)  Sd=2j+1  An.dQj+i)  ^  dlj  B2d-2j,  if  0  <  j  < 

.(2j  +  l)(2/),  if  j  =  m 


777- 


Five  initial  rows  of  triangle  generated  by  j  >  0,  0  <  j  <  777  are 
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(3.29) 


m  =  0 
m  =  1 
m  =  2 
m  =  3 
m  =  4 
m  =  5 


30 


-120 


-1386 


-14 


660 


140 


630 


2772 


Figure  8.  Triangle  generated  by  Amj,  j  >  0,  0  <  j  <  m,  sequences  A302971  (numerators  of 

Arn  j)  and  A304042  (denominators  of  Am,j). 

Note  that  starting  from  row  m  >  11  the  terms  of  Triangle  (|3.29)  consist  fractional  numbers,  for 
example,  Anq  =  800361655623,  6.  One  can  find  complete  list  of  the  numerators  and  denominators 
of  Am>j  in  OEIS  under  the  identifiers  A302971  and  A304042,  respectively,  see  [T7],[l8].  To  verify 
the  terms  that  definition  (3.28)  produces  one  should  refer  to  Mathematica  codc0  Hereby,  let  be 
theorem 


Theorem  3.30.  For  every  non-negative  integers  n  >  0  and  m  >  0  holds 


nv  =  < 


Y)  Dm(n,k)=  Y  Y  Amjkj(n-k)j , 


for  V  =  2m  +  1 


1 <k<n  l<k<nO<j<m 

Y  h Dm(n,k)=  Y  n  D  Amtjkj (n  -  k)j ,  for  V  =  2m, 

l<fc<n  0  <j<m 


1  <k<n 


where  Dm(n,k )  is  defined  by  (3.17). 

One  can  verify  results  concerning  above  theorem  (3.30)  via  Mathematica  cod^j  Therefore, 


theorem  (3.30)  answers  to  the  question  (2.23)  positively,  since  for  every  m  >  0  exists  a  triangle, 


generated  by  Dm(n,k ),  n  >  0,  0  <  k  <  n,  such  that  odd  power  n2m+1  can  be  reached  as  sum  of 
n-th  row  of  corresponding  triangle  over  I<fc<nor0<fc<w-1.  Sequences  A287326,  A300656, 
A300785  are  partial  cases  of  theorem  (3.30)  for  m  =  1,2,3,  respectively. 

3.1.  Properties  of  Dm(n,k )  and  Amj.  Here  we  show  a  few  properties  of  definition  Dm(n,k), 


some  of  them  correlates  with  properties  of  partial  case  D\(n,  k)  in  2.15 

(1)  Sum  of  Amj  over  j  in  range  0  <  j  <m  gives 


0<j<m 


Am..i  —  2 


_  r)2ra+l 


- 1, 


where  Amj  is  defined  by  (3.28). 


(2)  Similarly  to  property  (2.21 )  of  particular  case  D\(n ,  k),  items  of  {Dm(n,  k)}^=0,  m  >  0,  n  > 
0  is  symmetric,  i.e 

Dm(n,  k)  =  Dmin,  n  -  k), 

for  all  A:  :  0  <  k  <  n. 

(3)  By  property  ([2])  for  every  integer  n  >  0,  m  >  0  immediately  follows 

n2m+l  _  ^  Dm(n,k)  =  J2  Dm(ni  k) 


n 


(4)  For  every  m  >  0  the  Ar 

(5)  For  each  m  >  0 


1 <k<n  0<k<n—l 

are  terms  of  the  sequence  A002457, 


Am,  0  —  1 


‘ def_2_12.txt  -  Mathematica  code,  implementation  of  definition  (3.281,  [25  . 


'-expression_2_l.txt,  -  Mathematica  code,  implementation  of  theorem  (3.30) 
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(6)  Assume  that  n  <  0,  then  theorem  (3.30)|can  be  applied  as 

for  V  =  2m  +  1 


V 

n  = 


E  Dm(n,k), 

l</c<|n| 

I]  T^Dm(\n\,k),  for  V  =  2m 

<  l<^<|n| 


Property  3.31.  (Linear  Recurrence  of  Dm(n,k).)  For  every  integer  n  >  0  in  D\(n,k),  D2[n,k), 
D${n,k)  hold  the  recurrent  relations 

Di{n  +  l,k)  =  2D\(n,  k)  —  D\{n  —  1,  k) 

D2(n  +  2,k)  =  3D2(n  +  1,  k)  —  3D2(n,  k)  +  D2(n  —  1,  k) 

D3(n  +  3,k)  =  4T>3(n  +  2,  k)  -  6D3(n  +  1,  k)  +  4D3(n,  k)  -  D3(n  -  1,  k) 


Review  the  coefficient  Dm(n,k),  which  defined  by  m-order  polynomial.  R’s  well  known  fact  that 
high  order  finite  difference  A m+kPm(x)  =  0  for  every  x,  where  Pm{x)  is  m  order  polynomial  and 
k  >  0,  [8].  Recall  Binomial  Transform  of  sequence  an.  D.  E.  Knuth  [35]  has  introduced  the  binomial 
transform  by 

(3.32)  an  =  (?)  (-l)fc«fe, 

fc=o  ^  ' 


In  particular,  an  =  A nan,  therefore,  for  Dm(n,k )  we  have 

^t  +  1 


(3.33) 


Vf  >  m  :  Dm(n,  k)t  =  ^ 

i>o 


(-i  y 


3 


Dm(n  +  t  -  j,k) 


=  0, 


hereby,  it  gives  us  right  to  represent  recursively  every  value  of  Dm(n  +  r,k),  0  <  r  <  t  +  1  as 

(3.34)  (-1  )r+1Ct  +  1']Dm(n  +  t-r,k)  =  ^  (-1  )j+1  +  1'\  Dm(n +  t  -  j,k) 

'  r  '  iez>o/{r}  v  J  / 


In  particular, 

(3.35)  Dm(n  +  t,k)  =  ^(-l)i+1  (*  +  1)  Dm(n  +  t-j,  k) 

j> i  \  3  J 


□ 


Hereby,  let  be  theorem 


Theorem  3.36.  By  property  (3.31),  particularly,  from  expression  (3.35),  for  every  integer  t  >  m 
follows 


(3.37) 


(»  +  *)2m+1  =  EE(-l 

k> 1 j> 1 


^(T) 


Dm(n  +  t  -  j,k) 


Proof.  Direct  consequence  of  theorem  (3.30)  and  (3.35).  □ 

3.2.  Example  of  use  of  theorem  (3.30).  In  this  subsection  we  show  a  detailed  application  of 
theorem  (3.30)|  In  this  subsection  we  highlight  the  corresponding  terms  of  Ar 


m,j,  0  <  j  <  m  and 

T(n,k),  1  <  k  <  n  with  different  colors  to  be  more  easily  to  see  regularity.  Recall  existing  pattern, 
that  is  triangle  of  coefficients  Amj  defined  by  (3.28) 
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m  =  0 

1 

m  =  1 

1 

6 

m  =  2 

1 

0 

30 

m  =  3 

1 

-14 

0 

140 

m  =  4  1 

-120 

0 

0 

630 

m  =  5  1 

-1386 

660 

0 

0  2772 

Figure  9.  Triangle  generated  by  Amj, 

m  > 

0,  0  <  j  <  m. 

By  received  formula  Xu <k<n  ^j> o  ^m(«j  k)  =  n2m+1  each  line  of  above  triangle  being  multiplied 
by  r-7>0(n,  k)  :=  k^(n  —  k)i  and  summed  up  to  n  or  n  —  1  over  k  from  0  or  1,  respectively,  will 
result  odd  power  of  n2m+1,  depending  on  the  row  Amj,  0  <  j  <  mn  is  applied.  Consider  the  case 
n  =  3,  m  =  2,  we  introduce  triangle  built  using  T(n,  k )  :=  k(n  —  k),  n  >  1,  1  <  k  <  n  as  upper 
triangular  array,  let  be  1  <  n  <  5,  1  <  k  <  n 


n  =  1  n  =  2  n  =  3  n  =  4  n  =  5 


(3.38) 


3 

4 
3 
0 


4 

6 

6 

4 

0 


Figure  10.  Upper  right  triangle  generated  by  T(n,  k )  =  k{n  —  k),  n  >  1,  1  <  k  <  n,  sequence 

A094053,  E9|  in  OEIS. 


Example  3.39.  Consider  |theorem  (3.30)|  let  be  n  =  3  and  m  =  2,  then 

32'2+1  =  1  •  2°  +  0  •  21  +  30  •  22 

+  1  •  2°  +  0  •  21  +  30  •  22 

+  1  •  0°  +  0  •  01  +  30  •  02 

=  121  +  121  +  1  =  243 

Items  in  above  array  are  terms  of  the  third  row  of  triangle  A300656, 


Example  3.40.  Consider  theorem  (3.30)  let  be  n  =  3  and  m  =  3,  then 

32-3+i  =  i  .  2°  -  14  •  21  +  0  •  22  +  140  •  23 

+  1  •  2°  -  14  •  21  +  0  •  22  +  140  •  23 

+  1  •  0°  -  14  •  01  +  0  •  02  +  140  •  03 

=  1093  +  1093  +  1  =  2187 

Items  in  above  array  are  terms  of  the  third  row  of  triangle  A300785, 


Example  3.41.  Consider  theorem  (3.30)  let  be  n  =  4  and  m  =  3,  then 

42'3+1  =  1  ■  3°  -  14  •  31  +  0  •  32  +  140  •  33 

+  1  •  4°  -  14  •  41  +  0  •  42  +  140  •  43 

+  1  •  3°  -  14  •  31  +  0  •  32  +  140  •  33 

+  1  ■  0°  -  14  •  01  +  0  •  02  +  140  •  03 

=  3739  +  8905  +  3739  +  1  =  16384 
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Items  in  above  array  are  terms  of  the  forth  row  of  triangle  A300785.  We  can  perform  same  result 
for  42'3+1  in  terms  of  multiplication  of  certain  matrices, 

3°  31  32  33 
40  41  42  43 

3°  31  32  33 
0°  01  02  03 


^2-3+1  _ 


=  [1,  1,  1,  1]  * 


1 

-14 

0 

140 

Example  3.42.  Consider  theorem  (3.30)  let  be  n  =  5  and  m  =  3,  then 


-2-3+1  _ 


= 

1 

•  4°  — 

14- 

41  +  0  •  42  +  140  ■  43 

+ 

1 

•  6°  — 

14- 

61  +  0  •  62  +  140  •  63 

+ 

1 

•  6°  — 

14- 

61  +  0  •  62  +  140  •  63 

+ 

1 

•  4°  — 

14- 

41  +  0  •  42  +  140  ■  43 

+ 

1 

•  0°  — 

14- 

01  +  0  •  02  +  140  ■  03 

= 

8905  +  30157  +  30157  +  8905  +  1  =  78125 

Items  in  above  array  are  terms  of  the  fifth  row  of  triangle  A300785. 


Example  3.43.  Consider  theorem  (3.30)  let  be  n  =  5  and  m  =  4,  then 
52'4"1 


= 

1 

•  4°  — 

120 

•  41  +  0  •  42  +  0  •  43  +  630  •  44 

+ 

1 

•  6°  — 

120 

•  61  +  0  •  62  +  0  •  63  +  630  •  64 

+ 

1 

■  6°  — 

120 

•  61  +  0  •  62  +  0  •  63  +  630  •  64 

+ 

1 

■  4°  — 

120 

•  41  +  0  •  42  +  0  •  43  +  630  •  44 

+ 

1 

•  0°  — 

120 

•  01  +  0  •  02  +  0  •  03  +  630  •  04 

= 

160801  +  815761  +  815761  +  160801  +  1  =  1953125 

4.  Another  power  identity 


Review  the  Corollary  (2.2),  it  says  that  partial  case  of  theorem  (3.30)  returns  the  binomial  of 
the  form 


^2  Di(n,k) 


M<k<N 


f  U\{T,  0)n°  +  U\  (T.  1  )n\  if  M  =  1,  N  =  T 

|l7i(T-l,0)n0  +  l7i(T-l,l)n1,  if  M  =  0,  N  =  T-  1 

=  n3,  as  T  — >  n. 


Let  extend  this  idea,  as  we  have  complete  theorem  (3.30)  We  rewrite  theorem  (3.30)  for  odd 
powers  as 


^  '  Dm(n,  k) 

M<k<N 


J  Um(T,  0 )n°  +  •  •  •  +  Um(T,  m)nm,  if  M  =  1,  N  =  T 

\Um(T  —  1,  0)n°  +  •  •  •  +  Um{T  —  1,  m)nm,  if  M  =  0,  N  =  T-  1 

n2m+1,  as  T  =  n. 
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Above  expression  is  so-called  ’closed  form’  of  the  theorem  (3.30)  for  every  particular  n.  Below  we 
place  a  few  examples  of  polynomials  Xo<fc<m(— l)m kUm(T,  k )  •  nk ,  where  m  =  2, 


(4.1) 


T 

U2(T,  0)nu  +  U2(T,  l)nl  +  U2(T ,  2 )nz 

1 

31  —  60n  +  30n2 

2 

512  -  540n  +  150?z2 

3 

2943  -  2160n  +  420n2 

4 

10624  -  6000n  +  900n2 

5 

29375  -  13500n  +  1650n2 

6 

68256  -  26460n  +  2730n2 

7 

140287  -  47040n  +  4200n2 

8 

263168  -  77760n  +  6120n2 

9 

459999  -  121500n  +  8550n2 

10 

760000  -  181500n  +  11550n2 

Figure  11.  Table  of  polynomials  generated  by  Xa <k<T  D2{n,  k)  =  n5  as  T  — >  n,  see  definition 


fl3T7|)  for  D2(n,k). 

Below  we  show  a  plots  of  few  first  polynomials  Xa<fc<T  D2(n,  k)  for  T  =  1,2,3  and  compare  it 


with  corresponding  monomial  n5  by  the  theorem  (3.30) 


Figure  12.  Local  approximations  of  monomial  n 5  by  corresponding  polynomials  Xkakt  k) 


for  T  =  1,2,  3,  by  theorem  (3.30) 


We  can  see  that  monomial  n'5  can  be  easily  approximated  in  neighborhood  of  every  particular 
natural  n.  The  polynomials  from  figure  (|ll|)  can  be  generated  using  Mathematica  cod(f*|  To 
understand  the  nature  of  coefficients  Um(T.  k),  let  derive  them  directly  by  means  of  theorem  (3.30)[ 


(4.2) 


^2  ^  Am,jkJ(n  -ky  =  J2  Amfik°{n  —  k)°  4 - b  Am,mkm{n  -  k)r 

k=  1  j= 0  k= 1 


^2  Amflk°(n  -  k)°  +  ^2  Am,ik^ (n  -  k)1  H - b  ^  Am,mkm(n  -  k)m 

k=  1  k= 1  k=  1 

n  n  n 

Am,  0  ^2  k°(n  ~  fc)°  +  ^2^(11  -  k)1  ~\ - b  Am,m  ^2  km(n  -  kY 


k= 1 


k=  1 


k= 1 


=  n2m+1. 


4  Um(n,k)_coefficients2.txt.  -  Mathematica  code,  generates  the  polynomials  from  figure  (111,  [35].  If  necessary,  one 
could  re-define  m  >  0  within  the  code. 
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By  the  binomial  theorem, 


n—  1 


n—  1 


y~](n  —  k)3  kJ  =  EE(  ■) 

/c=0  /c=0  2—0 

We  rewrite  the  main  result  of  (|4.2|)  as 


(4.3) 


Amfl^2k°(n-  k)°  +  ylm,i^fc1(n  -  k)1  4 - -  &)T 

/c=l  fc=l  fc=l 

a,«E£(“)  n°-i(_l)^0+i  +  ^TOji  ee(T) 


fc=l 2=0 
n 


+  Am,2^Z  E  L>2  '( 


fc=l i=0 


fc=l  *=0 
n 


k= 1  *=0 


.^EE  ;  n3_*(— l)*fc3+* 


fc=l  2=0 


+  ^EE(T)  nm-i(_i)iA.m+i  =  n 


2m+l 


Rewrite  expression  (|4.3[)  in  extended  vie 

n 

(4.4)  An,0  ^  “  fc)°  +  X]  fcl(n  ~  H - ^  Am,m  ^2  ~  ky 


view 

n 


k= 1 


fc=l 


fc=l 


^4?rt,0  ^  ^ 


I  n°(— l)0/c0 


fc=i  L 

n  r 


+  Am,i  ^  (^H-l)0,1  +  f  n°{—l)lk2 


k= 1  L 


fc=l  L 

n  r 


+  (d^M^+nAt-ijv+QnVi)2^ 


+  ^m,3 


fc=l  L 


hn3(-l)“*3+  P)„2(_l)l*4+  Q„1(_1)V+  P')„»(-l)3fc» 


4“  Am}m  ^  ) 


fe=l  L 


™jnm(-l)°fcm  +  r  ^(-l)1^1  +  •  •  •  +  r  n°(-l)mfc2 


=  n2m+1. 
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Let  rewrite  expression  (|4.4[)  again  and  move  sigma  notation  under  the  brackets, 

n  n 

(4.5)  n2m+1  =  Amfi  E  k\n  -  k)°  +  ■  ■  •  +  Am,m  E  km(n  -  k)m 


k= 1 


fc=l 


=  -*4m,0 

+  Am2 

+  Am  3 


L  »"H)v 


Lk=l 
n 


+  A 


m,  1 


_k= i  v  '  fc=i  v  7 


Lfc=i 


k= i 


E  Ur2(“1)0fc2 + E  i j"1(-i)lfc3 + E  ( ;  w-i)2*4 


k= i 


Lfc=l 


fc=l 


E  o  ?z3(-1)0fc3 + E  ( *  M-i)1*4 + £  ( 2 Kt-1)2*5 + E  ( 3  A-1)3*6 


k= i 


fc=i 


+  -4„ 


E  nm(-l)°fcm  +  E  nm-1(-l)1fem+1  +  ---  +  ^  n°(-l)mfc2m 
Lfc=i  ^ 0 '  k= A 1 '  k,=\  w 

For  example,  consider  the  polynomial  (— l)m~fcC/m(n,  fc)nfc,  where  m  =  3.  We  derive  below  a 

set  of  coefficients  C/3 (n,  0),  C/3(n,l),  C/ 3(n,2),  1/3(71,3) 

3 

(4.6)  ^(-l)3-fcC/3(n,fe)nfc 

/c=0 


=  a 


3,0 


'O' 


E  hw 


Lfe=l 
n 


+  A; 


3,2 


+  -*4.3,1 


U=1 


EfMn^-D-P  +  ElLA-DV 


fe=l 


Lfc=i 


fc=i 


n2(-l)»t2  +  23  Wf-l)1*,3  +  23  ,  «°(-l)2P 


k= 1 


+  A; 


3,3 


=  n7. 


E  ( +  E  ( ?V2(-i)^4  +  E  uVm)2*5  +  E  uK(-i)3*6 


lk=\ 


k= 1 


fc=l 


k= 1 


The  coefficients  C/3 (n,  0),  C/3(n,  1),  C/3(n,  2),  C/3(n,  3)  in  expression  (4.6)  are  following 

n  /n\  n  /i\ 

(4-7)  C/3(n,0)  =  J43,oE(0)(-1)0fc0  +  243,iE(1j(-1)1/c2 


fe=i 

n 


k= 1 


fc=l 

n 


+  -43,2  E  (  9  )  (-1)2fc4  +  ^3,3  E  (  (— 1)3A)6, 


fc=l 


(4.8) 


C/3(n,l)  =  -43,i  E  Q  (-l)0*1  +  -43,2  E  (T)  (-!)^ 
fc=l  '  '  k= 1  '  7 


l/,3 


+  4.E  ,  (-DV. 


fe=l 
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(4.9) 


k= 1 


u3(n,2)  =  A3)2  (  j(-l)°fc2  +  A3i3y~]  (  )(-l)1fc' 


k= 1 


(4.10)  U3(n,  3)  =  ^3,3E(o)(“1)0/c3- 

fe=i  '  ' 

Let  rewrite  above  identities  in  compact  form  as, 

3  n 


(4.11) 


U3(n,r  =  0)  =  /)(-1) 


t>  o  £=1 

3  n 


\trilt 


(- 1)‘^ 


t>0  £=1 
3 


(4.12) 


C/3(n,r  =  l)  =  An,*  5])  (  )(-!)*  ^  1 


t>l  fcl 

3  n 


t- i 


t>i  f=i 
3  n 


(4.13) 


U3(n,r  =  2)  =  yi  Amit 

t>2  * 

3  n 

-  EE* 


t 


t> 2  £=1 

3  n 

t>2  £=1 
3 


t-2 

t 

t  -  2 


(4.14) 


U3(n,r  =  3)  = 


t>3  £=1 


(_l^-2£2*-2 
^ _ -Qt— 2  £2t—2 

\t-3n2t-3 


i=\ 


AaE  , v 


t  -  3 


Note  that  above  compact  forms  are  depend  on  r  =  0, 1,2,3,  as  it  can  be  seen  by  corresponding 
binomial  coefficients,  thus  we  have  replaced  k  by  £,  to  avoid  confusion  in  notations.  Therefore,  for 
every  0  <  r  <  m  the  Um(n,  r )  can  be  defined  next  way 


m  n 


(4.15) 


Um(n,  r )  =  ^2  Am,t  ^2 

t>r  £=  1 

m  n 

=  Ami 


t 

t  —  r 
t 

t  —  r 


(-l)t— r^t— r 

(- l)t-rft~r 


t>r  1=1 

Similarly,  the  coefficients  Um(n  —  1,  r)  derived  just  by  changing  the  iteration  limits  within  the  sum 
Ya=\  Am,t(2r)(-lY~r£2t~r  in  (4.15)  from  £  €  {l,n}  to  £  <E  {0,n  -  1} 

m  n—  1 

(4.16) 


um(n-  l.r)  e  E'4”>.>EL_J(-1)‘"r<2,"r 


t>r  1=0 

m  n—  1 

,,t 

t>r  1=0 


t 


^2  y2 


t  —  r 


t  —  r 


(—lyr£ 


t—r  [flt—r 
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In  above  equation  (4.16)  it  must  be  defined  =  1  for  all  i,  see  m-  Mathematica  implementations 
[^]of  expessions  (4.15)  and  (4.16)  are  available  in  [JTJ,  [42].  Therefore,  for  every  n  >  1,  and  m  >  0, 
we  have  identity 

(4.17)  n2m+1  =  < 

0  <r<m 

(4.18)  =  X]  ( 

0  <r<m 


-ly 

-ly 


Um(n,  r)  ■  nr 
" Um(n  -  1  ,r)  ■  nr 


Expressions  (4.17),  (4.18)  are  analogs  to  Faulhaber’s  odd  power  identity,  see 


,  p.  9.  Expressions  (4.17),  (4.18)  could  be  verified  via  Mathematica  codes 


property  (2.3)  and 
™  Here  we’d  like 


to  show  a  few  examples  of  polynomials  generated  by  X^o<fc<T-i  k),  see  definition  (3.17)  for 

D2(n,k) 


(4.19) 


T 

U2(T  -  1, 0)nu  +  U2(T  -  1,  IK  +  U2{T  -  1,  2 )nz 

1 

1 

2 

32  —  60n  +  30n2 

3 

513  -  540n  +  150n2 

4 

2944  -  2160n  +  420n2 

5 

10625  -  6000n  +  900n2 

6 

29376  -  13500n  +  1650?i2 

7 

68257  -  26460n  +  2730?i2 

8 

140288  -  47040n  +  4200n2 

9 

263169  -  77760n  +  6120n2 

10 

460000  -  121500?!  +  8550n2 

Figure  13.  Table  of  polynomials  generated  by  X^o<fc<T-i  ^(n,  k). 

Below  we  show  a  plots  of  few  first  polynomials  Xa<fc<T_i  ^(n,  k)  for  T  =  1,2,  3, 4  and  compare 


it  with  corresponding  monomial  n5  by  the  theorem  (3.30) 


Um(n,k)_coefficients_row_generating.txt.  -  Mathematica  code,  implementation  of  (4.151,  prints  ten  initial  rows  of 
coefficients  Um(n,  k) 


f  Um(n-l,k)_coefficients_row_generating.txt.  -  Mathematica  code,  implementation  of  (4.161,  prints  ten  initial  rows  of 
coefficients  Um(n  —  1  ,  k) 


'  Um(n,k)_odd_power_identity.txt  -  Mathematica  code,  implementation  of  ( 4.17 1,  [371 . 
Um(n,k)_odd_power_identity2.txt  -  Mathematica  code,  implementation  of  (4.181,  [38j 
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Y  Di(n,k) 


M<k<N 


Figure  14.  Local  approximations  of  monomial  n5  by  corresponding  polynomials 
Ei<KT-i^2(n,^)  for  T  =  1,2, 3, 4,  by  theorem  (3.30) 

The  polynomials  from  figure  (13)  can  be  generated  using  Mathematica  cod^J  Additionally,  the 
generalized  binomial  series  could  be  reached  by  means  of  identity 

f  U\  (T.  0 )n°  +  U\(T,  1 )n\  if  M  =  1,  N  =  T 

j[/i(T-l,0)n0  +  17i(T-  l,l)n\  if  M  =  0,  N  =  T  —  1 

=  n3,  as  T  — >  n. 

For  instance,  for  every  natural  n 

n3  =  Ui(T,  0)  •  n°  +  Ui{T,  1)  •  n1 
To  reach  higher  power  m,  let  just  multiply  this  identity  by  nm~3, 

nm  =  Ui(T,  0)  •  nm~2  +  U^T,  1)  ■  nm~3 
Let  rewrite  this  expression  regarding  to  itself  as  recursion, 

nm  =  Ui (T,  0)  •  [E/i  (T,  0)nm"4  +  Ux  (T,  l)nm~5}  +  Ui  (T,  1)  •  [E7i  (T,  0)nm"5  +  t/i  (T,  l)nm"6] 

=  L/i (T,  0)2  •  nm"4  -  2  ■  C/i  (T,  0)  •  t/i  (T,  1  )nm"5  +  C/i  (T,  l)2  •  nm“6 
Repeating  above  process  similarly  j  >  1  times,  we  have 


(4.20) 


n  = 


fc>0 


vfc 
xfc/j 


^(-!)fen  0)j"fc  •  U\(T,  l)fc  •  nm~2^~k 


k>  0 


2^(-!)fc  (^)  C/i(T  -  1,  0)J'-fc  -U\{T  —  1,  l)fc  •  nm~2j~k 


We  believe  we  can  perform  similarly  in  terms  of  multinomial  coefficients, 


(4.21)  n”  =  £ 

fci+fc2+---+fct=i 

=  £ 

fci+fc2+---+fct=i 


k\,  lv2,  ■  ■  ■  ,  kt 


k\i  &2j  •  •  •  j  kt 


f[(-l)eUt(T,£)h 


£=1 

t 


■  nrn-(t+1')i-k . 


l=i 


Note  that  we  must  set  T  =  n  in  above  identities.  Another  way  to  represent  odd-powered  monomial 
n2m+1,  m  =  0,1,2,...  is  to  define  a  polynomial  Um(n,k).  As  JJo<k<m(~ l)TO~fef/m(n,  k)  ■  nk  = 

—  1,  k)  ■  nk,  it  follows  that 


2  n 


2m+l 


Y  (-1  )m~kUm(n,k)-nk  +  Y  (-ir~kUm(n-l,k)-nk 

0<k<m  0  <k<m 

Y  (-1  r~k(Um(n,  k )  +  Um(n  -  1,  k))  •  nfc, 


0  <k<m 


Let  be  definition 

(4.22) 

Therefore, 

(4.23) 


Um{n,k )  d=  *  Um(n,  k)  +  Um(n  -  1,  k) 


n 


2m+l 


—  y  ^  kfmijli  k~) 


Tl 


0  <k<m 


"  Um(n,k)_coefficients.txt.  -  Mathematica  code,  generates  the  polynomials  from  figure  (12  I,  [TD] .  If  necessary,  one 
could  re-define  m  >  0  within  the  code. 
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Expression  (4.23)  generates  the  following  polynomials  given  m  =  2, 


(4.24) 


T 

U2{T,  0)nu  +  Z72(T,  1)?P  +  U2{T,  2)nz 

1 

16  -  30n  +  15 n2 

2 

272  -  300n  +  90n2 

3 

1728  -  1350n  +  285n2 

4 

6784  —  4080n  +  660n2 

5 

20000  -  9750n  +  1275n2 

6 

48816  -  19980n  +  2190n2 

7 

104272  -  36750n  +  3465n2 

8 

201728  -  62400n  +  5160n2 

9 

361584  -  99630n  +  7335n2 

10 

610000  -  151500n  +  10050n2 

Figure  15.  Table  of  polynomials  generated  by  XX <k<m^m(T,  k )  •  nk ,  where  T  =  1,2,  3, ..,  10. 
Polynomials  from  Figure  (15)  could  be  generated  using  Mathematica  codt|10l  Let  graphically  show 


an  approximation  of  monomial  n5  by  polynomials  XX <fc<m^2(T,  k)  ■  nk  for  T  =  1,  2,  3, 4. 


6784  -  4080  n  +  660  n2 
1728-  1350/1  +  285  n2 


272 -300n  + 90  n2 


16-30/1+  15/12 


Figure  16.  Local  approximations  of  monomial  n5  by  corresponding  polynomials 
XX<fc<m^2(T,  L)  •  nk  for  T  =  1,  2, 3, 4. 


Definition  (4.22)  of  coefficients  Um{n ,  fc)  and  identity  (4.23)  could  be  verified  via  Mathematica 
codes  m  respectively.  Similarly  to  (|4.20|),  (|4.2l|)  we  can  perform  a  binomial  and  multinomial 
representation  of  monomial  nm,  where  n,m  are  non-negative  integers  and  T  =  n, 


(4.25) 


nm  = 


k>  o 


^(-l)fc  J  )UAT,  0 y~k  ■  U\  (T,  l)fc  •  n 


m—2j—k 


k>  0 


J^(-l)fc(  ^  )^i(T-  l,0)J'“fc  -^i(T-  1, 1)*  .nm-2j~k 


^Combined_U_m(n,k)_coefficients_polynomials_gf.txt.  -  Mathematica  code,  generates  the  polynomials  from  Figure 

Mi, 


12 


JCombined_U_m(n,k)_coefficients_gf.txt.  -  Mathematica  code,  verifies  the  values  definition  (4.221  produces,  1441. 


Combined_U_m(n,k)_coefficients_odd_power_identity.txt.  -  Mathematica  code,  verifies  identity  (4.231,  [45] . 
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And 

(4.26) 


£ 

k1+k2+...+kt=j 


£ 

ki+k2+...+kt=j 


(  j 

\ki,k2,  ■ 


,  ^  T\(-l)eUt(T,i)kt  ■  nm-(t+1)j-k 
■'W  1=1 


)  IJ i-rfUtiT  -  1  J)k‘  ■ 
■  >  tJ  e=1 
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6.  Conclusion 


In  this  paper  particular  pattern,  that  is  binomial  distributed  triangle  A287326  in  OEIS,  which 
shows  perfect  cube  n  as  sum  of  row  terms  over  0  <  fc  <  n  —  1  or  1  <  fc  <  n  is  generalized,  in 
this  manuscript  are  found  and  discussed  the  polynomials  Dm(n,k )  and  Um(n,k),  such  that,  when 
being  summed  up  over  k  in  some  range  with  respect  to  m  and  n  returns  the  monomial  n2m+1. 
As  first  step,  we  discussed  analogs  of  A287326  for  powers  /  =  5,7,  sequences  A300656,  A300785, 
respectively,  then  we  derived  coefficients  Amj,  such  that  for  every  n  >  0  and  m  >  0  holds 

n2m+1  _  ^  Dm(n,  k), 

l<k<n 

where 


Dm(n,  k )  :=  Am,mkm(n  -  k)m  +  L(n  -  k)m  + 


+  Amfik{ 


\n-kf 


anc 


(3) 


Am,j  is 


defined  by  (3.28[).  Therefore,  question  question  (2.23)  is  answered  positively.  Section 


is  totally  dedicated  to  complete  and  extended  derivation  of  identity  ^ 


KKn  Dm(n,  k )  —  n 


_  „2m+l 


Properties  of  triangle  (2.14)  and  polynomial  Dm(n,k)  are  shown  in  properties  2.15  and  subsection 


3.1  respectively.  Relation  between  Faulhaber’s  sum  ^  nm  and  finite  differences  of  power  are  shown 
2.3  In  section  [4]  we  have  generalized  the  main  result  of  Corollary  (2.2)  for  all  odd  powers  of  the 


m 


form  2 m  +1,  m  =  0, 1,  2, ...  and  proven  an  identities 

n2m+l  =  ^  '  (-l)m~rU, 


£ 

0  <r<m 

£ 

0  <r<m 


(-i)r 


..  (n,  r)  ■  nr 
Um(n  -  1  ,r)  ■  nr 
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7.  Application  1.  Extended  tables  of  polynomials,  consisting  Um(n,k )  coefficients 

In  this  application  we  attach  an  extended  table,  consisting  of  polynomials,  generated  by  Yli<k<T  A 
and  So<fc<x— l  k)  for  various  t  and  m.  We  begin  from  cases  m  =  2,3,4  given  generating 

function  X)i<fc<T  Dm{n,  k)  and  continue  similarly  with  examples  for  X^o<fc<T-i  k).  The  fol¬ 

lowing  tables  could  be  generated  using  Mathematica  code  Um(n,k)_coefficients2.txt.  Here  we  begin 
to  show  our  tables  for  m  =  1,  2, 3, 4  and  T  =  1,2, ...,  40 


t 

Polynomial(n) 

1 

31  —  60n  +  30n^ 

2 

512  -  540n  +  150?r2 

3 

2943  -  2160n  +  420n2 

4 

10624  -  6000n  +  900n2 

5 

29375  -  13500n  +  1650n2 

6 

68256  -  26460n  +  2730n2 

7 

140287  -  47040n  +  4200n2 

8 

263168  -  77760n  +  6120n2 

9 

459999  -  121500n  +  8550n2 

10 

760000  -  181500n  +  11550n2 

11 

1199231  -  261360n  +  15180n2 

12 

1821312  -  365040n  +  19500n2 

13 

2678143  -  496860n  +  24570n2 

14 

3830624  -  661500n  +  30450n2 

15 

5349375  -  864000n  +  37200n2 

16 

7315456  -  1109760n  +  44880n2 

17 

9821087  -  1404540n  +  53550n2 

18 

12970368  -  1754460n  +  63270n2 

19 

16879999  -  2166000n  +  74100n2 

20 

21680000  -  2646000n  +  86100n2 

21 

27514431  -  3201660n  +  99330n2 

22 

34542112  -  3840540n  +  113850n2 

23 

42937343  -  4570560n  +  129720n2 

24 

52890624  -  5400000n  +  147000n2 

25 

64609375  -  6337500n  +  165750n2 

26 

78318656  -  7392060n  +  186030n2 

27 

94261887  -  8573040n  +  207900n2 

28 

112701568  -  9890160n  +  231420n2 

29 

133919999  -  11353500?r  +  256650n2 

30 

158220000  -  12973500?r  +  283650n2 

Figure  17.  Table  for  m  =  2,  generating  function:  Xm  <k<T  ^m(n,k)  over  T  =  1,2,. ..,30. 


t 

Polynomial  (n) 

1 

-125  +  406n  -  420n2  +  140n3 

2 

-9028  +  13818n  -  7140n2  +  1260n3 

(n,k) 
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3 

-110961  +  115836n  -  41160n2  +  5040n3 

4 

-684176  +  545860n  -  148680n2  +  14000n3 

5 

-2871325  +  1858290n  -  411180n2  +  31500n3 

6 

-9402660  +  5124126n  -  955500n2  +  61740n3 

7 

-25872833  +  12182968n  -  1963920n2  +  109760n3 

8 

-62572096  +  25945416n  -  3684240n2  +  181440n3 

9 

-136972701  +  50745870n  -  6439860n2  +  283500n3 

10 

-276971300  +  92745730n  -  10639860n2  +  423500?r3 

11 

-524988145  +  160386996n  -  16789080n2  +  609840n3 

12 

-943023888  +  264896268n  -  25498200n2  +  851760n3 

13 

-1618774781  +  420839146n  -  37493820n2  +  1159340n3 

14 

-2672907076  +  646725030n  -  53628540n2  +  1543500n3 

15 

-4267591425  +  965662320n  -  74891040n2  +  2016000n3 

16 

-6616398080  +  1406064016n  -  102416160n2  +  2589440n3 

17 

-9995653693  +  2002403718n  -  137494980n2  +  3277260n3 

18 

-14757360516  +  2796022026n  -  181584900n2  +  4093740n3 

19 

-21343778801  +  3835983340n  -  236319720n2  +  5054000n3 

20 

-30303773200  +  5179983060n  -  303519720n2  +  6174000n3 

21 

-42311023965  +  6895305186n  -  385201740n2  +  7470540n3 

22 

-58184203748  +  9059830318n  -  483589260n2  +  8961260n3 

23 

-78909220801  +  11763094056n  -  601122480n2  +  10664640n3 

24 

-105663629376  +  15107395800n  -  740468400n2  +  12600000n3 

25 

-139843308125  +  19208957950n  -  904530900n2  +  14787500n3 

26 

-183091507300  +  24199135506n  -  1096460820n2  +  17248140n3 

27 

-237330365553  +  30225676068n  -  1319666040n2  +  20003760ra3 

28 

-304794997136  +  37454030236n  -  1577821560n2  +  23077040n3 

29 

-388070250301  +  46068712410n  -  1874879580n2  +  26491500n3 

30 

-490130237700  +  56274711990n  -  2215079580n2  +  30271500n3 

31 

-614380739585  +  68298954976n  -  2602958400n2  +  34442240n3 

32 

-764704580608  +  82391815968n  -  3043360320n2  +  39029760n3 

33 

-945510081021  +  98828680566n  -  3541447140n2  +  44060940n3 

34 

-1161782683076  +  117911558170?r  -  4102708260n2  +  49563500n3 

35 

-1419139853425  +  139970745180n  -  4732970760n2  +  55566000n3 

36 

-1723889362320  +  165366538596n  -  5438409480?r2  +  62097840n3 

37 

-2083091040413  +  194491000018?r  -  6225557100n2  +  69189260n3 

38 

-2504622113956  +  227769770046n  -  7101314220?r2  +  76871340n3 

39 

-2997246219201  +  265663933080?r  -  8072959440n2  +  85176000n3 

40 

-3570686196800  +  308671932520n  -  9148159440n2  +  94136000n3 

Figure  18.  For  m  =  3,  generating  function:  Yli<k<T  Dm(n,  k)  over  T  =  1,2,. ..,40. 


t 

T 

2 

3 

4 

5 

6 

7 

8 


Polynomial(n) 

751  -  2640n  +  3780n2  -  2520n3  +  630n4 
162512  -  325440n  +  245700n2  -  83160n3  +  10710n4 
4297023  -  5837040n  +  3001320n2  -  695520?r3  +  61740n4 
45586624  -  47125200n  +  18484200n2  -  3276000n3  +  223020n4 
291683375  -  244000800n  +  77546700?r2  -  11151000n3  +  616770n4 
1349845776  -  949440240n  +  253906380n2  -  30746520n3  +  1433250n4 
4981676287  -  3024769440n  +  698619600n2  -  73100160n3  +  2945880n4 
15551330048  -  8309593440n  +  1689523920n2  -  155675520n3  +  5526360n 
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9 

42670773999  -  20362676400n  +  3698370900n2  -  304479000n3  +  9659790n4 

10 

105670786000  -  45562677600n  +  7478370900n2  -  556479000n3  +  15959790n4 

11 

240716895551  -  94670349840n  +  14174871480n2  -  962327520n3  +  25183620n4 

12 

511605381312  -  184966507440n  +  25461891000n2  -  1589384160n3  +  38247300n4 

13 

1025515755823  -  343092771840n  +  43707229020n2  -  2525042520n3  +  56240730n4 

14 

1955262884624  -  608734803600n  +  72168875100n2  -  3880359000n3  +  80442810n4 

15 

3569884005375  -  1039300430400n  +  115225437600n2  -  5793984000n3  +  112336560n4 

16 

6275713432576  -  1715757781440n  +  178643314080n2  -  8436395520n3  +  153624240n4 

17 

10670440655087  -  274981 1239440n  +  269883324900n2  -  12014435160n3  +  206242470n4 

18 

17613015856848  -  4292605722240ra  +  398449531620n2  -  16776146520n3  +  272377350n4 

19 

28312660615999  -  6545162506800?r  +  576282961800n2  -  23015916000?r3  +  354479580n4 

20 

44440660664000  -  9770762509200n  +  818202961800n2  -  31079916000n3  +  455279580n4 

21 

68269062114351  -  14309505635040n  +  1142398899180n2  -  41371850520n3  +  577802610n4 

22 

102840862500112  -  20595287515440n  +  1570974936300n2  -  54359003160n3  +  725383890n4 

23 

152176783290623  -  29175447644640n  +  2130550596720n2  -  70578587520n3  +  901683720n4 

24 

221524231290624  -  40733355636000n  +  2852919846000n2  -  90644400000n3  +  1110702600n4 

25 

317654602459375  -  56114215014000n  +  3775771408500n2  -  115253775000n3  +  1356796350?r4 

26 

449215653223376  -  76354376660640n  +  4943473041 780n2  -  145194842520?r3  +  1644691230?r4 

27 

627146261293887  -  102714466735440n  +  6407922490200n2  -  181354088160n3  +  1979499060n4 

28 

865161520339648  -  136716646589040n  +  8229467839320n2  -  224724215520n3  +  2366732340n4 

29 

1180316760605999  -  180186334891200n  +  10477899992700n2  -  27641231  lOOOn3  +  2812319370?r4 

30 

1593659760714000  -  235298734894800n  +  13233519992700?r2  -  33764831  lOOOn3  +  3322619370n4 

31 

2130981114417151  -  304630522458240n  +  16588283906880?r2  -  409793771520n3  +  3904437600?r4 

32 

2823673440038912  -  391217063149440n  +  20647028001600?r2  -  494350940160n3  +  4565040480n4 

33 

3709710869661423  -  498615539455440n  +  25528776924420n2  -  592972130520n3  +  5312170710n4 

34 

4834761029884624  -  630974381822400n  +  31368137616900n2  -  707469399000n3  +  6154062390?r4 

35 

6253442526125375  -  793109409951600n  +  3831 678 1679400n2  -  839824524000n3  +  7099456140n4 

36 

8030741767978176  -  990587103477840ra  +  46545018909480n2  -  992199287520n3  +  8157614220?r4 

37 

10243603824112687  -  1229815433857440n  +  56243464735500n2  -  1166946059160n3  +  9338335650n4 

38 

12982712871538448  -  1518142701993840n  +  67624804267020?r2  -  1366618682520n3  +  10651971330n4 

39 

16354478705823999  -  1863964838829600n  +  80925655683600?r2  -  1593983664000n3  +  12109439160n4 

40 

20483246706016000  -  2276841638834400?r  +  96408535683600n2  -  1852031664000n3  +  13722239160n4 

Figure  18.  For  m  =  4,  generating  function:  Xacfccr  k)  over  T  =  1,2,. ..,40. 
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